The multi-scale dust formation in substellar 
atmospheres 

Christiane Helling 1 ' 2,3 , Rupert Klein 2,4 ' 5 , and Erwin Sedlmayr 1 

^^ ' x Zentrum fiir Astronomie und Astrophysik, TU Berlin, Hardenbergstrafie 36, 

D- 10623 Berlin 

sedlmayrOastro . physik . tu-berlin . de 
<*"] ' 2 Konrad-Zuse- Zentrum fiir Informationstechnik Berlin, Takustrafie 7, D- 14195 

^ , Berlin 

rupert . kleinOzib . de 

3 Sterrewacht Leiden, P.O Box 9513, 2300 RA Leiden 
helling@strw . leidenuniv . nl 

4 Fachbereich Mathematik und Informatik, Freie Universitat Berlin, Takustrafie 7, 
D- 14195 Berlin 

5 Potsdam Institut fiir Klimafolgenforschung, Telegrafenberg A31, D-14473 
Potsdam 

.«—■>. ■ Summary. Substellar atmospheres are observed to be irregularly variable for which 

^—4- ' the formation of dust clouds is the most promising candidate explanation. The 

f^ , atmospheric gas is convectively unstable and, last but not least, colliding convective 

"*!> ■ cells are seen as cause for a turbulent fluid field. Since dust formation depends on 

the local properties of the fluid, turbulence influences the dust formation process 
and may even allow the dust formation in an initially dust-hostile gas. 
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© ' A regime-wise investigation of dust forming substellar atmospheric situations 

reveals that the largest scales are determined by the interplay between gravitational 
settling and convective replenishment which results in a dust-stratified atmosphere. 
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The regime of small scales is determined by the interaction of turbulent fluctuations. 
Resulting lane-like and curled dust distributions combine to larger and larger struc- 
tures. We compile necessary criteria for a subgrid model in the frame of large scale 
simulations as result of our study on small scale turbulence in dust forming gases. 



1 Introduction 

The astrophysical field of substellar objects - of which brown dwarfs shall be 
considered as an example - has gained considerable attention during the last 
few years since more and more brown dwarfs and extrasolar planets could 
be detected. The need for appropriate models has increased with increasing 
observational resolution power which, however, can only provide information 
about the largest scale structures. The major challenge has been the obser- 
vational evidence for the presence of dust (i.e. small solid or fluid particles) 
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in substellar atmospheres where energy is mainly transported by convection. 
These gas parcels collide after their ascend through the atmosphere during 
which they adiabatically expand. Such adiabatic collisions cause turbulence, 
and the turbulent kinetic energy is transfcred cascade-like from large to small 
scales. Due to its chemical nature, the dust formation process itself (like com- 
bustion) depends on local quantities like temperature, density and chemical 
composition of the gas, which can vary on much smaller spatial scale lengths 
than those accessible by observations. Photon interactions, molecular colli- 
sions and friction influence the dust formation process and proceed on scales 
comparable to or even less than one mean free path of the fluid, which belongs 
to the smallest relevant scale regime. The dust formation is hence initiated 
and controlled by small-scale processes, but produces consequences on the 
largest, observable scales. 

The immediate couplings between chemistry, hydrodynamics and thermo- 
dynamics due to the presence of dust can cause an amplification of initially 
small perturbations resulting in the formation of large-scale dust clouds. The 
formation of such clouds is one of the major candidates to explain the ob- 
served non-periodic variability [?, ?, ?, ?, ?, ?, ?] which is - from the physical 
point of view - very alike to what we know from weather-like variation in the 
Earth atmosphere. However, substellar object like brown dwarfs and known 
extraterrestrial planets show more extreme physical conditions (warmer and 
hotter or cooler and thinner) and the direct transfer of Earth's knowledge has 
to be considered with care. 

It was therefore the aim of this work to investigate the multi-scale prob- 
lem of astrophysical dust formation in brown dwarf stars which requires an 
appropriate description of the chemical and turbulent processes and their in- 
teractions. The challenge in modeling turbulence in reactive gas flows lies in 
an adequate description of all relevant scale regimes. This work follows the 
scale hierarchy from the smallest to the largest scales. At first, turbulent dust 
formation will be studied on small scales (Sect. 3) which can be computation- 
ally resolved. These investigations intend to provide the basis for a sub-grid 
closure model for the next (larger) scale regime where it can be applied and so 
on. This procedure requires a detailed, scale-dependent view on the problem 
[?], or in other words, needs to adopt different windows of perception [?]. 



2 The model of a dust forming gas flow 

Redefining the variables of the system of model equations like a — > a/a le { 
transforms the equations in their scaled analogues wherein all quantities are 
dimensionless and can be compared by number. The non-dimensional char- 
acteristic numbers describe the importance of the source terms and general 
features can already be derived by investigating the dimensionless numbers 
for physical meaningful reference values. 
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Complex A: The specification of the source terms in the equations of momen- 
tum and energy conservation results in the following dimensionless system of 
equations which describe the fluid field in a substellar atmosphere which is 
influenced by the stellar radiation field [?] : 



(p) t +V-(pv) = (1) 

Sr(pv) t + V-(pvov) = —y S VP--pj S pg (2) 

{pe) t + V • (v[pe + P]) = Mi «(T* E - T 4 ), (3) 

with the caloric equation of state 

( pv 2 1 \ P 

pe = lM * [^ + — pgyj + — - (4) 

with g = {0, g, 0}. Radiative heating/cooling is treated by an relaxation ansatz 
(r.h.s. Eq. 3). For the characteristic numbers M, Fr, Sr, and Rd\, which relate 
certain physical processes, see Table 2. 

Complex B: The dust formation process is considered as a two step process. 
At first, seed particles form out of the gas phase (nucleation) which provide 
the first surfaces. Subsequent growth by surface reactions results in the for- 
mation of (chemically) macroscopic /mi-sized particles. Nucleation, growth, 
evaporation, drift and element depletion/enrichment are physical and chemi- 
cal processes which occur simultaneously in an atmospheric gas flow and may 
be strongly coupled. Following the classical work of [?], partial differential 
equations which describe the evolution of the dust component by means of 
moments of its size distribution function, f(V), were derived in [?] explicitly 
allowing for u gas ^ fdust ■ The conservation form of the dust model equations 
allows for a fast numerical solution in the frame of extensive model simula- 
tions: 



d_ 
dt 



oc oo 

{pLfi+V^pLj) = Jj2 R kV j/3 dV - yjf(V)V^ 3 v dT (V)dV,(5) 
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with v dr (V) the grain size dependen equilibrium drift velocity (for more details 
please consult [?]). The j-th moment of the dust size distribution function 
Lj [cm J /g] is defined by 

oo 

pLj(x,t) = I f(V,x,t)V^ 3 dV . (6) 
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The source term Aj expresses the effects of nucleation and surface chemical 
reactions on the dust moments. Compared to the classical moment equations 
[?], Bj is an additional, advective term in the new dust moment equations 
which comprises the effects caused by a size-dependent drift motion of the 
grains. The meaning of the dust variables is summarizes in Table 1 . 



Table 1. Definitions, meanings, and units of the dust and chemical quantities 
(CE = chemical equilibrium). 



Quantity Unit 



Meaning 



Variables 



t^gas 


[cms -1 ] 


hydrodynamic gas velocity 


^dust 


[cms -1 ] 


hydrodynamic dust velocity 


f(V) 


[cm" ] 


grain size distribution function 


V 


[cm 3 ] 


volume of the dust particle 


Jh 


K 1 ] 


surface chemical reaction rates 


h 


[cm j g 1 ] 


dust moments (j 6 N) 


pLo = rid 


[cm -3 ] 


number of dust particles 


Y 4-7T Lq \ / 


[cm] 


mean radius of dust particles 


v^xf = <^> 


[cm 2 ] 


mean dust surface 


% = (V) 


[cm 3 ] 


mean dust volume 


e x 


H 


element abundance relative to hydrogen 


n T 


[cm 3 ] 


number density of gas species r 


n ■ p 


[cm- 3 ] 


total hydrogen density 


"w 1.427 amu 


^rcl.x 


[cms -1 ] 


thermal relative velocity 
of the gas species x 



Constants 



Vifi 
Vi.r 



stoichiometric ratios of homogeneous nucleation 
stoichiometric ratios of surface reaction r 



The motion of a dust grain with a velocity iMust is determined by an equi- 
librium between the force of gravity, -F gra v, and the frictional force, -Ff r i C , 
(equilibrium drift => v dr (V); see discussion in [?]). Depending on the particle 
size and the density of the surrounding fluid, the character of the hydro- 
dynamic situation changes which also influences the growth process of the 
particles. The dimensionless dust moment equations for nucleation, growth, 
evaporation, and equilibrium drift write, e. g. for the case of a subsonic free 
molecular How, are 



[Sr] (pL j ) t + V(v gas pLj) = [Sr • Da7 c • Sej] J{V t ) 



(7) 



S r ' Da d,lKn 



Xx&pLj-i 



/7r7\V2 Sr-M-Dr 
V32/ 



6k„V 



l j+i 

C T 



e r 
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with j — 0, 1, 2 . . . the order of the dust moments. The first source term on the 
r.h.s. describes the seed formation (nucleation) , the second the growth process, 
and the third describes the transport of already existing dust particles by drift. 
In those cases, where dust and gas are positionally coupled and drift effects 
are negligible, Eq. (7) simplifies to 

[Sr] (pL j ) t + V(v gas pLj) = [St. Dag™ • Sej] J(V t ) + [Sr- Da| r ] | x not pL^ x . 

(8) 

For the definition of the characteristic numbers in square brackets see Table 2. 
The element depletion is taken into account by evaluating the consump- 
tion of each involved element i with relative abundance to hydrogen ej by 
nucleation and growth 

R 
[Sr] (pe x ) t + V ■ (v pe x ) = - [Sr- El] ^«f Da£ uc N L J* 

i-i 

+ i/& t Da| r a r n Xjr w rc i jX pL 2 ). (9) 

e x is the element abundance of the chemical element x (e.g. Ti, Si, O) in mass 
fractions, which is consumed by nucleation (first term, r.h.s.) and growth 
(second term, r.h.s.) with the stoichiometric factor v x ^ T due to the reaction 
r. One equation needs to be solved for each element involved in the dust 
formation process. 

Characteristic behavior 

Adopting typical values for a brown dwarf atmosphere the following general 
characteristics of a dust forming fluid can be drawn (for definitions see Ta- 
ble 2): 

High Reynolds numbers, Re ss 10 7 . . . 10 9 indicate that the viscosity in the 
brown dwarf atmosphere is too small to damp hydrodynamical perturba- 
tions and a turbulent hydrodynamic field can be expected. 
Assuming a typical turbulence velocity of the order of one tenth of the 
velocity of sound leads to a Mach number of M sw O(0. 1) for the whole 
fluid. Fluctuations on smaller scales can have M « 0(1O -2 ). 
The Froude number Fr = 0(1O -3 . . . 10 _1 ) shows that the gravity only 
gains considerable influence on the hydrodynamics for scale regimes Z re f ^ 
H p . Drift terms are important in the macroscopic regime and may be ne- 
glected in the small scale regimes. An analysis of the combined character- 
istic drift number shows that the drift term in the dust moment equations 
is mainly influenced by the gravity and the bulk density of the grains. 



Table 2. Characteristic numbers and reference values of the scaled model equations in Sect. 2. 



S' 



o 
x 






cr 

bjj 
= 
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Characteristic numbers 



Mach number 

Froude number 

Strouhal number 

hydrodyn. Knudsen number 

Knudscn number 

Drift number 

1. Radiation number 

Sedlmayr number (j £ No) 

Damkohler no. of nucleation 

Damkohler no. of growth (lKn^ 

Element consumption number 



M = 

Fr= ■ 
Sr = 

HD 



~cf 1 

of Bret 

'rcf 

rofrcf 

'rcf 

2<a>rcf 

Kn= J^- 

2 Wrcf 

Dr = ^^£ 



Kn 



rof ' P rof 
j 



Rdi = 4 Krcf aT r 4 ef ■ 

So 



Da 1 



Da „UC = trcf^cf 
fl Prof-to.rcf 

_ <rofXrof,lKn 
d.lKn _ 7E 

El = 



(W^f) 1 ^ 

Piof-C-O.rcf"; 



II i.l'i-l 



Combined characteristic number 



combined drift number (lKn) 



7T7\i/2 Sr-M-Dr 



32/ 



Kn HO - Fr 



Independent reference values 
(fundamental dimensions) 



density 


Pref 


temperature 


T re { 


velocity 


Vrel 


length 


Iret 


ravitational acceleration 


ffref 


nucleation rate 


^rof/^<H>,rof 


mean particle radius 


(«) rcf 


element abundance 


£ref 



Dependent reference values 



thermal pressure 

total hydrogen density 

hydrodynamic time 

th dust moment (= n^/p) 

growth velocity 

mean free path 

opacity 

molecular number density 



[g/cm 
[K 
[cm/s 
[en 
[cm/s 



cm 

2 

1 



18 

[cm 



D Picf kT icf 






[dyn/cm 2 


W<H>,rof = Pici/J2 e i m i 






[cm" 3 


i AV Si0 n ral Si0 , t s 

-Lo.rof fin/a)* { > 
Fret 3 \"/ re f 

XlKn 






[s 

[1/g 

[cm/s 


4cf 






[cm 


/ abs.eas r \ 
K re f = K ro f(o- rot B ,1^3, rof) 






[cm- 1 


n re f = n rofj sio ~ ej^<H>,rcf 


c 


) 


[cm" 3 



lKn — large Knudscn number {Kn ^> 1) (*) — approximation (compare Eq. 5) ^**^ _ ^ Q ^ c determined from chemical equilibrium calculations 



Parameters: a - Stcfan-Boltzmann constant 
pd.rcf - bulk density [gem - ], 77 



L s cm' K 4 



ai - hypothetical monomer radius [cm], /AVsiO - monomer volume of SiC>2 [cm 



element mass [g] , 0"p e / - gas absorption cross section 
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The characteristic number for the radiative heating / cooling, Rdi = 
4K re fcrT r g f • -pst- gives the ratio of the radiative and the thermal energy 
content of a fluid. The scaling of the system influences Rdi by the refer- 
ence time i re f which increases with increasing spatial scales. 
An analysis of the characteristic numbers in front of the source terms of 
the dust moment equations reveals a clear hierarchy of the processes of 
dust formation: nuclcation — > growth — ► drift. 

The governing equations of the model problem are those of an inviscid, 
compressible fluid which are coupled to stiff dust moment equations and an 
almost singular radiative energy relaxation if dust is present. The coupling 
between the processes becomes stronger with increasing time scales. 



If additionally v gas — > 0, i. c. the dust-forming system reaches the static case, 
and t — > oo, the source terms in Eqs. (7) must balance each other. In the 
S > 1 case (S - supersaturation ratio), this means that the gain of dust by 
nucleation and growth must be balanced by the loss of dust by rain-out. In 
the S < 1 case, just the opposite is true, i. e. the loss by evaporation must 
be balanced by the gain of dust particles raining in from above. Both control 
mechanisms (in the static limit) result in an efficient transport of condcnsiblc 
elements from the cool upper layers into the warm inner layers, which cannot 
last forever. 

If the brown dwarf's atmosphere is truly static for a long time, there is 
no other than the trivial solution for Eqs. (7) where the gas is saturated 
(S=l) and dust- free (Lj = 0). 



3 Dust formation on small scales 
3.1 Numerical approach 

The fully time-dependent solution of the model equations (Eqs. 1 - 3, 8, 9) 
has been obtained in the small scale regimes by applying a multi-dimensional 
hydro code [?] which has been extended in order to treat the complex of dust 
formation and elemental conservation. 

Initial conditions: 

The initial conditions have been chosen as homogeneous, static, adiabatic, 
and dust free, i. e. po — 1, Po = 1, Uq = 0, Lq = (=> Lj = 0) in order to 
represent a (semi-)static, dust-hostile part of the substellar atmosphere. This 
allows us to study the influence of our variable boundaries on the evolution of 
the dust complex without a possible intersection with the initial conditions. 
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Boundary conditions and turbulence driving: 

The Cartesian grid is divided in the cells of the test volume (inside) and the 
ghost cells which surround the test volume (outside). The state of each ghost 
cell is prescribed by our adiabatic model of driven turbulence for each time. 
The hydro code solves the model equation in each cell (test volume + ghost 
cells) and the prescribed fluctuations in the ghost cells are transported into 
the test volume by the nature of the HD equations. The numerical boundary 
occurs between the ghost cells and the initially homogeneous test volume and 
are determined by the solution of the Riemann problem. Material can flow 
into the test volume and can leave the test volume. The solution of the model 
problem is considered inside the test volume. 

Stiff coupling of dust and radiative heating / cooling 

Dust formation occurs on much shorter time scales than the hydrodynamic 
processes. Approaching regimes of larger and larger scales makes this problem 
more and more crucial. Therefore, the dust moment and element conservation 
equations (Complex B) arc solved applying an ODE solver in the framework 
of the operator splitting method assuming T, p =const during ODE solution. 
In [?] we have used the CVODE solver (Cohen & Hindmarsh 2000; LLNL) 
which turned out to be insufficient for the mesoscopic scale regime. CVODE 
failed to solve our model equations after the dust had reached its steady state. 

The LIMEX solver: 

The solution of the equilibrium situation of the dust complex is essential for 
our investigation since it describes the stationary case of Complex B when 
no further dust formation takes place. The reason may be that all available 
gaseous material has been consumed and the supersaturation rate 5 = 1 
or the thermodynamic conditions do not allow the formation of dust. The 
first case involves an asymptotic approach of the gaseous number density of 
5=1 which often is difficult to be solved by an ODE solver due to the 
choice of too large time steps. However, the asymptotic behavior is influenced 
by the temperature evolution of the gas/dust mixture which in our model is 
influenced by radiative heating/cooling. Since the radiative heating/cooling 
(Qrad = Rdn(T^ E — T 4 ) heating/cooling rate) depends on the absorption 
coefficient n of the gas/dust mixture which strongly changes if dust forms. 
Consequently, the radiative heating/cooling rate is strongly coupled to the 
dust complex which in turn depends sensitively on the local temperature which 
is influenced by the radiative heating/cooling. It was therefore necessary to 
include also the radiative heating/cooling source term in the separate ODE 
treatment for which we adopted the LIMEX DAE solver. 

LIMEX [?] is a solver for linearly implicit systems of differential algebraic 
equations. It is an extrapolation method based on a linearly implicit Euler 
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discretization and is equipped with a sophisticated order and stepsize control 
[?]. In contrast to the widely used multistep methods, e. g. OVODE, only linear 
systems of equations and no non-linear systems have to be solved internally. 
Various methods for linear system solution are incorporated, e. g. full and band 
mode, general sparse direct mode and iterative solution with preconditioning. 
The method has shown to be very efficient and robust in several fields of 
challenging applications in numerical [?, ?] and astrophysical science [?]. 

3.2 Microscopic regime 

The smallest scale involved into the structure formation of a substellar atmo- 
sphere is the monomer volume, AV, by which a chemical reaction r increases 
the volume of a grain. For radiative transfer and drift effects, the grain size 
becomes important. On small hydrodynamic scales, drift effects are neglig- 
gibaly small. However, the aim is to investigate the role of turbulence in the 
formation of dust clouds and the possibly related variability of brown dwarfs. 
Therefore, the interesting hydrodynamic phenomena in the microscopic scale 
regime arc acoustic waves. 



dust ^r^HM *^^. radiative 

formation (CT r\ cooling 




dominance of dust opacity 



Fig. 1. Unstable feedback loop in dust forming systems: Dust, once formed, increases 
the total opacity which intensifies radiative cooling. The temperature decreases and 
enters the dust formation window. Therefore, more dust is formed which increases 
the opacity further. An instability establishes where the dust improves the conditions 
for its own formation f?l. 



Carrying out numerical simulations of the whole time-dependent system of 
model equations (Eqs. 1 - 3, 8, 9), turbulence is seen as acoustic waves being 
by-products of colliding convective cells [?] which undergo an adiabatic in- 
crease of size during their upward travel through the substellar atmosphere. A 
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feedback loop establishes in a dust-forming system in which interacting expan- 
sion waves play a key role (Fig. 1): The interaction of small-scale disturbances 
of the fluid field can cause a local and temporarily limited temperature de- 
crease low enough to initiate dust nuclcation. These seed particles grow until 
they reach a size where the dust opacity is large enough to accelerate radia- 
tive cooling which causes the temperature to decrease again below a nucleation 
threshold (T s ). Dust nuclcation is henceforth re-initiated which results in a 
further intensified radiative cooling. The nucleation rate and consequently 
also the amount of dust particles increase further. This run-away process is 
stopped if either the radiative equilibrium temperature of the gas is reached 
or all condensible material is consumed. Meanwhile, the seed particles have 
grown to macroscopic sizes. 

The result is a highly variable dust distribution in space and time on the 
microscopic scales investigated (e. g. Al f=a 10 2 cm, At w 0.5 s) due to the 
occurrence of singular nucleation events. It is immediately clear, that such 
scales are not easily resolvable by observations, but these are the scales on 
which structure formation is likely to be initiated. However, such small scale 
pattern might move and successively initiate dust formation at various sites in 
the atmosphere. A larger and larger cloud structure may thereby form which 
is, hence, determined by a non-local hydrodynamic coupling whereas at each 
site the above outlined local feedback loop will act (compare Fig. 3). 

3.3 Mesoscopic regime 

Following Kolmogoroff's idea, the mesoscopic scales shall be considered as 
those where energy is only transfered through the turbulence cascade, but is 
not injected or dissipated. Convective elements are generated as long as the 
Schwarzschild criteria is fulfilled in the substellar atmosphere. The inertia of 
mass will drive the convective elements beyond Schwarzschild's boundary, a 
phenomenon usually named as overshooting, which extends this zone further. 
The mass elements interact, collide, and transfer part of their energy into 
small scale structures, thereby producing a whole spectrum of them being 
viewed as waves or - turbulence elements. These waves propagate and enter 
even atmospheric regions which are not influenced by convective motions any 
more. 

The question arises how the dust formation process takes place in such a 
stochastically fluctuating thcrmo- and hydrodynamic environment, and there- 
fore which influence such small scale disturbances may have on the whole at- 
mospheric structure and the formation of large scale pattern possibly respon- 
sible for observed variabilities. A model for driven turbulence was proposed 
[?, ?] which simulates a constantly occurring, small-scale energy input as- 
sumed to originate from a convectively active zone. Following Reynolds scale 
separation ansatz a background field ao(x,t) is disturbed by a fluctuation 
5a(x, t) of some variable a(x, t) such that 
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a(x, t) — ao(x, t) + 6a(x, t). 



(10) 



a(x, t) e {u(x, t),P(x, t), S(x, £)} (u(x, t) - velocity, P - pressure, S - entropy). 
The velocity fluctuation 5u(x,t) follows the Kolmogoroff spectrum in fc-space 
in which a whole range of wavcnumbers is excited (k min = 2ir/l max , k max = 
2tt/(3Ax), / max = 5 x 10 4 cm - maximum scale considered, Ax = 10 2 cm - 
spatial grid resolution). The model for driven turbulence relies on a superpo- 
sition ansatz of different wave modes and therefore allows again to carry out 
direct simulation of the time and space evolution of the dust complex, now in 
a stochastically excited medium typical for the clswise dust-hostile part of a 
brown dwarf atmosphere. 
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Fig. 2. Fluctuating dust component in substellar atmospheres. The space- means 
(a(t)) x (solid/red) with the apparent standard deviations 0n x -i(*) (dotted/dashed) 
as function of time for T ra i — 2100K, p ro f = 3.16 10 _4 gcm~ 3 , v IC { — cs/10. ({a(t)) x + 
on x _i(£) - dotted/blue; (a(t)) x — 0n x -i(*) - dashed/green). 
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Nucleation events and nucleation fronts 

An initially hydrodynamically homogeneous and dust-free, 500m-sizcd gas 
parcel has been considered in a deep layer of a brown dwarf atmosphere ini- 
tially too hot for nucleation. The gas is disturbed by superimposed waves 
entering through its boundaries. If these disturbances carry already a temper- 
ature below the nucleation threshold (T s - nucleation threshold temperature) , 
a nucleation front will develop: The nucleation peak moves inwards together 
with the wave and initiates the feedback loop already known from the micro- 
scopic regime by leaving behind first dust seeds. If the entering temperature 
disturbances is not enough to cross T s , interaction with some expansion wave 
coming from another direction at some time and some site will take place 
thereby causing singular nucleation events to occur. 

A nucleation front tends to homogeneously fill the gas parcel with dust 
in ID situations while nucleation events cause a more heterogeneous dust 
distribution due to their short life time. In more than ID, also nucleation 
fronts cause a heterogeneous dust distribution because their 'parent' waves 
will interact and influence the conditions for dust formation. 

Large variations in the dust quantities occur when the formation process 
starts to influence the hydro- and thermodynamics of the gas parcel. If all 
available material has been consumed, the mean dust quantities reach an 
almost constant values (Fig. 2) at a certain place in the atmosphere which 
is characterized by T re f. During this active time (At rs 3s), the variation 
in the number of dust particles is O(10 5 cm -3 ) and in the mean particle 
size 0(lO 2 /im). The variation of the density is considerable (see standard 
deviations in Fig. 2) during the time of temperature decrease which assures 
the establishment of the pressure equilibrium. 

A convectively ascending, initially dust free gas element can be excited 
to form dust by waves running through it. A cloud can therefore be fully 
condensed at much higher temperatures than classically expected, i. e. in 
an undisturbed case. 

Formation of large scale structures in 2D 

2D simulations with the smallest eddy size being A™ n = 5 m, the largest are of 
the size of the gas box A^ x = 500 m, reveal a very intermittent distribution of 
the dust in space during the time of efficient dust formation (Fig. 3). The very 
inhomogeneous appearance of the dust complex is a result of nucleation fronts 
and nucleation events comparable to the ID results. The nucleation is now 
triggered by the interaction of eddies coming from different directions. Large 
amounts of dust are formed and appear to be present in lane-like structures 
(large logn<i; dark/red areas in Fig. 3). The lanes are shaped by the constantly 
inward traveling waves. Our simulations show that some of the small scale 
structure merge thereby supporting the formation of lanes and later on even 
larger structures. 
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Fig. 3. Spatial appearance of the number of dust particles (lognd [cm ]; false 
color background) and the vorticity (V x u; black and grey contour lines) of the 2D 
velocity field for t = 0.8s of a simulation with T re f = 2100 K, p rc f = 3.16 10 _4 gcm -3 , 

Vref = Cs- 



Dust is also present in curled structure which indicates the formation of 
vortices due to the 2D waves driving turbulence. As the time proceeds in the 
2D simulation, vortices develop orthogonal to the velocity field which show a 
higher vorticity (V x u(x, t)) than the majority of the background fluid field. 
Comparing the distribution of dust particles (false color background in Fig. 3) 
with the vorticity (contour lines in Fig. 3), shows that the vortices with high 
vorticity preferentially occur in dust free regions or regions with only little 
amounts of dust present. These vortices can efficiently drag the dust into 
regions with still material available for further condensation shaping thereby 
larger and larger structures. 



4 Dust formation on large scales 



The atmospheric extension (& 100 km - several pressure scale heights H p ) 
proposes a natural upper limit for the macroscopic scale regime. On this 
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scale, the interplay between convective overshooting and gravitational settling 
(drift) is a major mechanism which influence the dust structure of the whole 
atmosphere. Since the dust is strongly coupled to the hydro- and thermody- 
namics of the atmosphere, gravitational settling and convective overshooting 
indirectly alter the atmospheric density and temperature structure which are 
important for the spectral appearance of the objects. The large scale struc- 
tures are those which are usually accessible by observations. Therefore, the 
understanding and modeling of processes directly connected with macroscopic 
scale motions are of particular interest. 

4.1 The model of a quasi-static dust layer 

A theoretical consistent description of dust formation, gravitational settling, 
and element depletion (Eqs. 7, 9) was needed in order to understand the 
formation and the structure of dust cloud layers which is urgently needed for 
explaining the spectral appearance of Brown Dwarfs [?] . A first application to 
the quasi-static case of a brown dwarf atmosphere calculation was carried out 
by providing a model of a quasi-static cloud layer, i.e. v g3uS = and -gf- = 
in Eqs. (1) - (3), (7), (9) [?]. However, if v gas = Eqs. (7) have the trivial 
solution Lj = (see Sect. 2). The physical interpretation of this solution is that 
dust grains have once formed in the sufficiently cool layers, have consumed 
all available condensible elements up to the saturation level, and have finally 
left the model volume by gravitational settling. Consequently, a truly static 
atmosphere must be dust-free which - in this generality - contradicts the 
observations. 

This picture changes, however, if we take into account a mixing of the 
atmosphere caused by the convection. This mixing leads to an ongoing re- 
plenishment of the gas with fresh, uncondensed matter from the deep interior 
of the brown dwarf, which can counterbalance the loss of condensible elements 
via the formation and gravitational settling of dust grains. Thus, a convective 
mixing can maintain a stationary situation: Seed particles nucleate in metal- 
rich, i. c. supersaturated upwinds. Dust particles grow on top of these nuclei 
by the accretion of molecules and rain out as soon as they have reached a 
certain size. The sinking grains will finally reach deeper atmospheric layers 
which are hot enough to cause their thermal evaporation, which completes the 
life cycle of a dust grain in a brown dwarf atmosphere. We therefore enlarge 
the static case of our model by a simple description of the convective mixing 
such that the stationary dust moment equations write 

-T-(^) - 7^ i-^ + V^ + i^pL^) (11) 

dz\ C T J £lKn V r mix 3 / 

with algebraic auxiliary conditions 

n(H)(£t ~ frt) = v^NtJ* + ¥M^pL 2 f^^ r n r vfa r ( 1 - i-) . (12) 



r=l 



I J I 
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According to this approach, the gas/dust mixture in the atmosphere is con- 
tinuously replaced by dust-free gas of solar abundances on a depth-dependent 
mixing time-scale r m i x (z), which can be adapted to a convection model (see 
discussion in [?, ?]). q is the actual abundance of element i in the gas phase 
and e° its solar value [?] . By fitting the mass exchange frequency of 3D dy- 
namic simulations of surface convection by [?] with an exponential we can 
apply a rough description of T m i x (z). 

4.2 Numerical approach 

Equation (11) is a system of ordinary differential equations of first order which 
can be solved by standard numerical methods. The differential equations are 
integrated inward by means of the variable transformation z' = z max — z us- 
ing the Radau 5 - solver for stiff ordinary differential equations [?] ) . For test 
purposes, the equation of hydrostatic equilibrium dP glis /dz = —pg is solved 
in addition to the three moment equations (11), where the actual gas pres- 
sure P = ^2 i riikT results from the chemical equilibrium calculations. The 
integration is stopped as soon as one of the dust moments becomes negative, 
indicating that the dust has completely evaporated. 

4.3 Structure of a quasi-static cloud layer 

The vertical structure of the cloud layer results from a competition between 
the four relevant physical processes: mixing, nucleation, growth/evaporation 
and drift. Following the cloud structures inward (from the left to the right in 
Fig. 4) roughly five different regions can be distinguished, marked by the Ro- 
man digits, which are characterized by different leading processes concerning 
the dust component. 

0. Dust-poor depleted GAS: High above the convection zone, the mixing 
time-scale is large and the elemental replenishment of the gas is too slow to al- 
low for considerable amounts of dust to be present in the atmosphere. The few 
particles forming here are very small (a) < 10 -25 /im and have drift velocities 
(u dr ) rj 1 mm/s to 1 cm/s which causes these atmospheric layers to become 
dust-poor. The gas phase is strongly depleted in condensible elements. The 
Ti abundance in the gas phase exi is reduced by a factor between 10 5 and 
10 7 from its solar value. However, phase equilibrium (S = 1) is not achieved 
because even the very small disturbance of the atmosphere by mixing is suffi- 
cient to produce a solution which differs significantly from the trivial solution 
in the truly static case. 

1. Region of efficient nucleation: Nucleation takes place mainly in the 
upper parts of the cloud layer, where the temperatures are sufficiently low and 
the elemental replenishment by mixing is sufficiently effective. Although the 
gas is strongly depleted in heavy elements in these layers, it is nevertheless 
highly supersaturated (S'>1000) such that homogeneous nucleation can take 
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Fig. 4. Calculated Ti02 cloud structures for models with log g — 5 and T c g — 
1800 K. Note the scaling of the partly logarithmic and partly linear p-axis. 1 st panel: 
prescribed gas temperature T (solid), according to Tsuji (2002), and mixing time- 
scale T m i x (dashed) . 2 nd panel: nucleation rate J* (solid) and supersaturation ratio 
StiC^ (dashed) . 3 rd panel: Ti abundance in the dust phase e d " st (solid) and in 
the gas phase eii (dashed-dotted). The solar value e^i i s additionally indicated by 
a thin straight line. The dashed line shows the growth velocity XiKn- 4 th panel: 
mean particle size (a) = iJ/3/(47r) Li/Lo (solid) and mean drift velocity (vdr) = 
\ZTigpd{a} I (2pcr) (dashed, see Eq. (66) in Paper II). 5 th panel: molecular particle 
densities of TiO (solid) and Ti02 (dashed). For comparison, the hypothetical total 
titanium nuclei density for solar abundances n^ H )£Ti is depicted by a thin solid line. 
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place efficiently. Since very many seed particles are produced in this way, the 
dust grains remain small (a) < 0.01 /an and have small mean drift velocities 
(u dr ) w 0.1 mm/s ... 1 mm/s, which are even smaller than in region because 
of the higher gas densities. 

II. Dust growth region: With the inward increasing temperature, the su- 
persaturation ratio S decreases exponentially which leads to a drastic decrease 
of the nucleation rate J*. Consequently, nucleation becomes unimportant at 
some point, i.e. the in-situ formation of dust particles becomes inefficient in 
region II. Here, the condcnsiblc elements mixed up by convective overshoot 
are mainly consumed by the growth of already existing particles, which have 
formed in region I and have drifted into region II. The gas is still strongly 
supersaturated S^$> 1, indicating that the growth process remains incomplete, 
i. e. the condensible elements provided by the mixing are not exhaustively 
consumed by growth. The dust component in region II is characterized by 
an almost constant degree of condensation (oc e^" st ex L 3 sa const), while the 
mean particle size (a) and the mean drift velocity (udr) increase inward. Con- 
sequently, the total number of dust particles per mass (Lo) and their total 
surface per mass ((xL 2 ) decrease. 

III. Drift dominated region: With the decreasing total surface of the 
dust particles (0C-L2), the consumption of condcnsiblc elements from the gas 
phase via dust growth becomes less effective. At the same time, due to the 
increase of the mean particle size (a), the drift velocities increase. When the 
dust particles have reached a certain critical size, (a) C r~15/im to 50 /im, the 
drift becomes more important than the growth, and the qualitative behavior 
of the dust component changes. This happens at the borderline between re- 
gion II and III, which we denote by rain edge. Although the gas is still highly 
supersaturated 5~1 ... 10, the in-situ formation of dust grains is ineffective as 
in region II. The depletion of the gas phase vanishes in region III, i. c. the gas 
abundances approach close-to-solar values. The grains reach their maximum 
radii at the lower boundary of region III (the cloud base): (a) w 30 /xm to 
400 /im at maximum drift velocities. 

IV. Evaporating grains: The gravitationally settling dust particles finally 
cross the cloud base and sink into the undcrsaturated gas situated below, 
where S < 1 . Here, the dust grains raining in from above evaporate thermally. 
The evaporation of these dust particles, however, docs not take place instan- 
taneously, but produces a spatially extended evaporation region IV with a 
thickness of about 1 km. With decreasing altitude, the particles get smaller 
d(a)/dz < and slower d(vdr)/dz < 0. Consequently, their residence times 
increase, and a run-away process sets in which finally produces a very steep 
decrease of the degree of condensation, terminated by the point where even 
the biggest particles have evaporated completely. We note that region IV, in 
particular, cannot be understood by stability arguments, but requires a kinetic 
treatment of the dust complex. 
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5 Conclusions 

The dust formation in brown dwarfs atmospheres is a multi-scale problem 
where different physical mechanisms characterize the different scale regimes: 

a) In the microscopic scale regime, single acoustic waves interact, thereby 
initiate dust formation. 

b) The superposition of a whole spectrum of turbulence elements determine 
the mesoscopic scale regime resulting in a highly fluctuating fluid held. 

c) The interplay of gravitational settling of dust grains and the convective 
up-mixing of uncondensed material governs the macroscopic scale regime. 
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Fig. 5. Regimes of turbulent dust formation. Ts: nucleation threshold temperature 
(supersaturation S <C 1 required), T BU ^: sublimation temperature. 



The main scenario envisioned is a convectively ascending fluid element in 
a brown dwarf atmosphere, which is excited by turbulent motions and just 
reaches sufficiently low temperatures for condensation. The dust formation 
in a turbulent gas is found to be strongly influenced by the existence of a 
nucleation threshold temperature Tg. The local temperature T must at least 
temporarily decrease below this threshold in order to provide the necessary 
supersaturation for nucleation, e.g. by eddy interaction. 

Depending on the relation between the local mean temperature T and Ts, 
three different regimes can be distinguished (see l.h.s. of Fig. 5): (i) the deter- 
ministic regime (T < Ts) where dust forms anyway, (ii) the stochastic regime 
(T > Ts) where T < Ts can only be achieved locally and temporarily by tur- 
bulent temperature fluctuations, and (iii) a regime where dust formation is 
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impossible. The size of the stochastic regime depends on the available turbu- 
lent energy. This picture of turbulent dust formation is quite different from the 
usually applied thermodynamical picture (r.h.s. in Fig. 5) where dust is sim- 
ply assumed to be present whenever T < T su b, where T su b is the sublimation 
temperature of a considered dust material. 

After initiation, the dust condensation process is completed by a phase of 
active particle growth until the condcnsiblc elements are consumed, thereby 
preserving the dust particle number density for long times. However, radiative 
cooling (as follow-up effect) is found to have an important influence on the 
subsequent dust formation, if the dust opacity reaches a certain critical value. 
This cooling leads to a decrease of T(t) which may re-initiate the nucleation. 
This results in a runaway process (unstable feedback loop) until radiative and 
phase equilibrium is achieved. Depending on the difference between the initial 
mean temperature T(t = 0) and the radiative equilibrium temperature Ire, 
a considerable local temperature decrease and density increase occurs. Since 
the turbulent initiation of the dust formation process is time-dependent and 
spatially inhomogeneous, considerable spatial variations of all physical quan- 
tities (hydro-, thermodynamics, dust) occur during the short time interval of 
active dust formation (typically a few seconds after initiation), which actu- 
ally creates new turbulence. Thus, small turbulent perturbations have large 
effects in dust forming systems. Our 2D simulations show that the dust ap- 
pears in lane- like and curled structures. Small scale dust structures merge and 
form larger structures. Vortices appear to be present preferentially in regions 
without or with only little dust. Non of these structures would occur without 
turbulent excitation. 

From our work on small-scale turbulence simulations of dust- hostile regions 
in substellar atmospheres, we compile necessary criteria for a subgrid model 
of a dust forming, turbulent system. 
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Criteria on a subgrid model of a turbulent, dust-forming system: 

a) The subgrid model should describe the transition stochastic — ► deter- 
ministic regime in dust forming turbulent fluids. 

b) The dust formation process (nucleation + growth) is restricted to a short 
time interval since the dust formation time scales are much smaller than 
the large-scale hydrodynamic time scales. 

This involves that: 

• The nucleation does occur only locally and event-like in very narrow 
time slots. 

• The growth process continues as long as condensible material is avail- 
able. 

• The condensation process freezes in and the inhomogeneous dust 
properties are preserved. 

• Almost constant characteristic dust properties result in the mean-long 
term behavior. 

c) The feedback loop with its fast radiative cooling should govern the tran- 
sition from an almost adiabatic to an isothermal behavior of the dust/gas 
mixture. 

The largest, observable scale of a brown dwarf atmosphere has been inves- 
tigated by attacking the drift problem. A consistent theoretical description 
was derived for dust formation and destruction, gravitational settling, and cl- 
ement depletion including the effect of convective overshoot. It was therewith 
possible for the first time to overcome the widely used concept of ad hoc dust 
presents in substellar atmospheres. Test calculations have shown that the dust 
will appear stratified in the atmosphere causing a corresponding depletion of 
the gas phase. 
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